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This  work  presents  the  development  of  a  dynamic  SOFC-GT  hybrid  system  model  applied  to  a  long-haul 
freight  locomotive  in  operation.  Given  the  expectations  of  the  rail  industry,  the  model  is  used  to  develop 
a  preliminary  analysis  of  the  proposed  system’s  operational  capability  on  conventional  diesel  fuel  as  well 
as  natural  gas  and  hydrogen  as  potential  fuels  in  the  future.  It  is  found  that  operation  of  the  system  on  all 
three  of  these  fuels  is  feasible  with  favorable  efficiencies  and  reasonable  dynamic  response.  The  use  of 
diesel  fuel  reformate  in  the  SOFC  presents  a  challenge  to  the  electrochemistry,  especially  as  it  relates  to 
control  and  optimization  of  the  fuel  utilization  in  the  anode  compartment.  This  is  found  to  arise  from  the 
large  amount  of  carbon  monoxide  in  diesel  reformate  that  is  fed  to  the  fuel  cell,  limiting  the  maximum 
fuel  utilization  possible.  This  presents  an  opportunity  for  further  investigations  into  carbon  monoxide 
electrochemical  oxidation  and/or  system  integration  studies  where  the  efficiency  of  the  fuel  reformer  can 
be  balanced  against  the  needs  of  the  SOFC. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Long-haul  locomotive  transportation  is  one  of  the  most  preva¬ 
lent  methods  of  shipping  goods  in  the  United  States.  An  exploratory 
committee  formed  in  2002  at  the  Center  for  Transportation 
Research  at  Argonne  National  Laboratory,  under  the  guidance  of  the 
U.S.  Department  of  Energy,  reported  that  in  1997, 4  billion  gallons  of 
diesel  fuel  were  consumed  by  locomotives  in  the  United  States.  This 
amount  represents  10%  of  all  diesel  used  for  transportation  and 
2.3%  of  all  transportation  fuels.  For  the  rail  operators,  this  repre¬ 
sents  an  annual  cost  of  $2  billion,  contributing  7%  of  their  total 
operating  expenses  [1].  Today,  the  situation  is  not  much  different, 
as  the  consumption  of  petroleum  in  freight  rail  represented  2.1%  of 
the  total  national  transportation  fuel  use  in  2005,  thereby 
consuming  571.4  trillion  BTUs  of  fuel  [2].  Thus,  fuel  use  for  long- 
haul  locomotives  is  continually  of  economic  interest  to  both 
government  and  industry,  and  represents  a  field  with  much 
potential  for  increasing  energy  independence. 

In  addition  to  these  motivating  factors,  the  use  of  diesel  fuel  in 
railway  applications  presents  a  significant  environmental  concern. 
Although  the  nationally  averaged  emission  signature  of  locomotives 
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nationwide  appears  low  (in  2001,  locomotive-sourced  NO*  was  only 
5%  of  the  national  total),  locomotives  are  responsible  for  much  of  the 
emissions,  especially  diesel  PM,  in  the  areas  where  they  are 
stationary  [1].  At  classification  railyards,  locomotives  are  respon¬ 
sible  for  96%  of  PM  emissions;  at  intermodal  railyards,  the  contri¬ 
bution  is  39%  [3].  The  general  public  living  near  these  stationary  rail 
operations  is  most  affected  by  these  emissions,  and  improvements 
in  locomotive  emission  signatures  can  enhance  the  air  quality  for 
this  subset  of  the  population.  Moreover,  in  California,  rail  operations 
tend  to  be  centered  in  areas  that  have  difficulty  achieving  air  quality 
standards  due  to  high  source  concentration  and  stagnating  meteo¬ 
rological  conditions;  thus,  the  additional  burden  is  placed  on  those 
living  in  areas  already  known  to  have  poor  air  quality.  Current  esti¬ 
mates  attribute  2980  premature  deaths,  5100  cases  of  acute  bron¬ 
chitis,  62,000  cases  of  asthma  and  other  respiratory  problems,  and 
830  heart-related  hospital  admissions  per  year  to  goods  movement 
in  the  Southern  California  Air  Basin  [4].  Similar  to  passenger  vehicles 
and  diesel  freight  trucks,  locomotives  have  thus  been  assigned 
a  progressively  limiting  tier-based  schedule  of  emissions  standards. 
This  Tier  system  requires  substantial  ongoing  reductions  in  the 
emission  of  NOx,  CO,  hydrocarbons,  smoke,  and  particulates  [5]. 

Satisfying  the  ever-evolving  Tier  standards  may  require  inves¬ 
tigation  of  technology  beyond  the  conventional  combustion 
reciprocating  engine.  Early  investment  in  emerging  technologies, 
such  as  Solid  Oxide  Fuel  Cells  (SOFC),  may  provide  more  options 
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and  enable  development  of  systems  with  greater  environmental 
benefit.  The  fuel  cell  has  already  been  identified  as  having  the 
capability  to  replace  the  entire  diesel  and  electric  power  system  for 
a  locomotive,  simply  based  on  economic  and  performance  concerns 
[1].  Moreover,  the  proven  high  efficiency  and  low  emissions 
signature  [6,7]  of  SOFC-based  systems,  for  both  criteria  pollutants 
and  greenhouse  gases  like  CO2,  present  a  compelling  opportunity 
for  vast  improvement  in  the  overall  economic,  environmental,  and 
health  impacts  of  freight  railroading. 

In  spite  of  the  fact  that  fuel  cells  in  general  have  been  considered 
in  the  literature  for  the  freight  locomotive  application,  very  little 
work  has  been  presented  for  the  use  of  the  solid  oxide  fuel  cell  in 
particular.  Most  theoretical  studies  and  the  only  major  demon¬ 
stration  project  to-date  have  focused  on  utilizing  low- temperature 
proton  exchange  membrane  (PEM)  fuel  cells  [8-14].  Systems 
assuming  on-board  power,  off-board  power  with  electrified  line, 
and  hybridization  with  large  auxiliary  battery  storage  have  been 
investigated  [9,11-14].  The  last  of  these  is  the  subject  of 
a  substantial  demonstration  project  in  the  United  States  as  part  of 
a  collaborative  effort  between  the  Burlington  Northern  Santa  Fe 
Railway  System  and  Vehicle  Projects,  LLC  [15-17]. 

Although  the  PEM  systems  have  shown  promise,  their  operation 
requires  a  hydrogen  infrastructure,  a  constraint  that  may  be  too 
limiting  at  the  moment  for  widespread  adoption  and  which  the 
SOFC  can  avoid  thanks  to  its  fuel-flexible  operation.  Additionally, 
the  SOFC  can  be  hybridized  with  a  gas  turbine  (termed  an  SOFC- 
GT),  allowing  for  a  number  of  synergistic  benefits.  Such  a  system 
has  been  shown  to  potentially  operate  at  efficiencies  above  60%, 
produce  virtually  zero  levels  of  criteria  pollutants,  and  operate  with 
low  noise  levels  due  to  the  small  number  of  moving  parts.  The 
ability  to  accomplish  all  of  these  and  operate  on  the  currently  used 
diesel  fuel  is  a  capability  that  is  unique  to  the  combination  of  high- 
temperature  fuel  cell  with  gas  turbine. 

Since  the  introduction  of  the  SOFC-GT  concept,  there  has  been 
significant  advancement  in  understanding  prototypical  systems 
through  work  simulating  the  performance  of  various  designs  and 
exploring  detailed  considerations  such  as  system  integration, 
control  optimization,  and  drivers  and  challenges  in  the  expected 
dynamic  response  [18-27].  In  addition,  the  available  literature  is 
rich  with  a  variety  of  study  methodologies  and  system  design 
options  to  develop  in-depth  understanding  of  key  system-wide 
dynamic  interactions  and  component  limitations  as  well  as  iden¬ 
tification  of  preferred  system  designs.  Finally,  the  SOFC-GT  research 
community  has  also  been  able  to  demonstrate  the  potential  of 
these  systems  to  be  flexible  in  application,  through  analyses  of 
systems  applied  to  such  varied  goals  as  grid  stabilization,  Integrated 
Coal  Gasification-Fuel  Cell  plants,  and  Distributed  Energy  Resources 
[28-30]. 

However,  these  prior  research  works  are  all  based  on  the 
designs  of  SOFC-GT  systems  for  stationary  power.  With  its  high 
fuel-to-electric  conversion  efficiency,  the  SOFC-GT  is  an  ideal 
candidate  for  stationary  power  applications,  especially  large-scale 
and  centralized  power.  This  feature  is  also  attractive  in  the  loco¬ 
motive  application,  since  the  wheels  are  run  by  electric  traction 
motors.  Still,  the  combination  of  the  SOFC-GT  with  this  mobile 
platform  has  yet  to  make  an  appearance  in  the  literature.  Mobile 
SOFC-GT  systems  have  been  proposed  by  Winkler  and  Lorenz,  but 
for  the  smaller  passenger  vehicle  platform  [31  ].  There  is  also  some 
available  literature  regarding  the  application  of  SOFC-GT  systems 
for  marine  applications  [32,33].  As  of  the  time  of  this  publication, 
there  does  not  seem  to  be  any  investigation  of  the  potential  for  the 
SOFC-GT  hybrid  system  for  the  locomotive  application. 

Significant  capital  investment  in  new  fueling  infrastructure 
required  by  PEM  solutions  can  be  avoided  by  taking  advantage  of 
the  SOFC-GT  fuel  flexibility.  This  study  has  adopted  a  scenario  in 


which  diesel  fuel  is  the  near-term  choice,  followed  by  natural  gas  as 
a  bridging  fuel,  and  finally  hydrogen  as  an  end-goal  fuel  in  the 
future.  It  is  assumed  that  the  two  fossil  fuels  are  pre-reformed  off 
board  from  the  locomotive.  Thus,  this  work  analyzes  the  perfor¬ 
mance  of  the  system  when  operated  on  a  diesel  autothermal 
reformate,  a  natural  gas  steam  methane  reformate,  and  humidified 
hydrogen.  Possibilities  for  on-board  reformation  will  be  the  subject 
of  future  work. 

To  develop  the  model,  this  work  builds  upon  the  foundation  laid 
out  by  the  previous  studies  and  aims  to  provide  insight  into  the 
unique  issues  of  the  SOFC-GT  in  a  large  mobile  platform.  The  major 
issues  include  the  need  for  fuel  processing  to  accommodate  diesel 
fuel  use  and  the  prediction  of  SOFC  degradation  due  to  coke 
development  in  the  anode.  For  this  reason,  the  current  investiga¬ 
tion  develops  a  model  of  an  SOFC-GT  system  that  operates  on 
hydrogen  or  fuel  reformates,  but  is  flexible  enough  to  later  incor¬ 
porate  the  requirements  of  the  reformer  and  the  undesired  but 
prevalent  effects  of  coking  reactions.  Because  these  chemical 
processes  are  often  complex  and  include  reactions  on  vastly 
different  timescales,  the  model  was  developed  with  flexible  control 
over  numerical  algorithms  specifically  designed  for  the  solution  of 
stiff  ODEs  in  concert  with  simulation  of  the  remaining  pertinent 
physics  in  the  fuel  cell. 

Therefore,  a  novel  model  framework  was  developed  and 
implemented  in  FORTRAN,  which  leaves  open  the  possibility  of 
incorporation  of  various  numerical  schemes  for  reaction  chemistry 
necessary  for  investigating  the  dependency  of  system  performance 
on  fuel  choice.  In  addition,  since  future  versions  of  the  model  are 
projected  to  incorporate  a  voltage  loss  mechanism  based  on  the 
simulation  of  coking  deposition  within  the  fuel  cell,  the  method  of 
convergence  for  the  current-voltage  relationship  has  been  designed 
to  accommodate  for  more  than  just  the  traditional  loss  mecha¬ 
nisms,  which  often  allows  for  simplifications  in  the  SOFC  model 
development.  This  paper  presents  the  development  of  this  model 
and  insights  into  steady-state  fuel  dependence. 

2.  Modeling  methodology 

2.1.  Global  solution  method 

The  FORTRAN  model  developed  in  this  work  consisted  of  a  set  of 
modular  subsystems  based  on  the  analysis  of  each  major  compo¬ 
nent  in  the  system.  Each  subsystem  required  iteration  to  conver¬ 
gence  within  itself,  and  in  some  cases  required  repeated 
convergence  iteration  among  subsystems.  The  schematic  of  the 
system  modeled  in  this  work  is  as  shown  in  Fig.  1.  It  is  a  simplified 
model  of  the  extensive  system  expected  to  be  required  of  an 
SOFC-GT  powering  a  locomotive  and  includes  the  major  compo¬ 
nents  and  features  required  in  this  application.  Namely,  the  system 
includes  direct  linking  between  the  SOFC  and  gas  turbine  units, 
includes  an  auxiliary  combustor  before  the  turbine,  includes  an  air 
preheater  supported  by  system  waste  heat,  and  includes  cathode 
recycle  to  aide  SOFC  thermal  management.  Some  balance  of  plant 
components  are  not  accounted  for  in  this  model  (such  as  fuel 
pumps,  a  fuel  preheater,  and  any  necessary  fuel  processing  equip¬ 
ment);  thus,  the  performance  found  in  this  investigation  will  be 
slightly  improved  compared  to  a  more  detailed  system  model.  It  is 
expected  that  this  simplification  still  provides  a  reasonable 
expectation  of  the  system’s  capabilities. 

The  overall  system  model  is  divided  into  four  major  subsystems: 
the  SOFC  subsystem,  the  gas  turbine  subsystem,  the  air  preheater 
subsystem,  and  the  controller.  The  first  three  of  these  are  closely 
tied  via  fluid  streams  passing  through  all  of  these  subsystems.  Thus, 
convergence  within  each  subsystem  must  be  matched  by  global 
convergence  of  values  shared  between  subsystems.  In  order  to 
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Fig.  1.  Base  SOFC-GT  system  simulated  and  analyzed  in  this  work. 


achieve  this,  the  gas  turbine  subsystem  was  chosen  for  the  overall 
convergence  check.  Given  the  shaft  speed  and  pressure  ratio 
utilized  for  interpolation  on  the  compressor  performance  map,  the 
turbine  model  calculated  an  expected  mass  flowrate  based  on 
similar  interpolation.  Agreement  between  this  flowrate  and  the 
flowrate  leaving  the  combustor  determined  overall  system 
convergence.  Within  the  air  preheater  and  the  SOFC,  internal 
convergence  criteria  determined  the  gas  and  component  states  for 
a  given  global  pair  of  gas  turbine  shaft  speed  and  pressure  ratio.  If 
global  system  convergence  was  not  obtained,  the  pressure  ratio  was 
adjusted  and  all  subsystems  reevaluated.  In  order  to  close  the 
solution  at  the  air  preheater,  the  current  timestep  air  from  the 
compressor  and  previous  timestep  exhaust  from  the  turbine  were 
utilized.  Accuracy  was  maintained  due  to  the  small  timesteps 
utilized  in  this  investigation. 


2.2.  SOFC  subsystem  model 

The  solid  oxide  fuel  cell  subsystem  is  based  on  a  quasi  two- 
dimensional  dynamic  model,  with  spatial  resolution  along  the  flow 
path  of  the  anode  and  cathode  gases  for  a  single  fuel  cell.  It  is 
assumed  in  the  model  that  the  anode  and  cathode  flows  are 
arranged  in  a  co-flow  orientation.  An  entire  fuel  cell  stack  is 
assumed  to  be  comprised  of  multiple  cells,  each  of  which  has 
multiple  channels  with  the  same  conditions  as  those  solved  for  by 
the  fuel  cell  module.  This  simplification  is  illustrated  in  Fig.  2, 
depicting  the  reduction  from  a  stack  of  cells  to  a  representative 
central  cell,  and  further  to  a  representative  central  channel  pair 
(anode  and  cathode)  with  associated  solid  structures  (fuel  cell 
tri-layer  and  current  collector).  Finally,  in  order  to  achieve  spatial 
resolution  in  the  model,  this  representative  geometry  is  then  sub¬ 
divided  along  its  flow  direction  into  nodes.  In  this  work,  each 
channel  was  simulated  with  10  nodes;  the  total  area  of  the  simu¬ 
lated  square  cell  was  0.25  m2.  Each  node  contained  four  control 
volumes:  anode  gas,  cathode  gas,  tri-layer,  and  current  collector.  In 
addition,  given  the  small  length  and  timescales  of  the  simulation, 
each  node  was  assumed  to  behave  similar  to  an  individual  perfectly 
stirred  reactor.  This  allows  each  node  in  the  channel  to  be  solved 
sequentially,  with  the  state  of  a  given  node  acting  as  the  inlet  of  the 
next  node. 

The  SOFC  nodes  are  each  solved  in  the  same  manner,  through 
the  completion  of  five  major  processes:  the  thermodynamics, 
conservation  of  momentum,  reaction  kinetics,  electrochemistry 
and  conservation  of  mass,  and  heat  transfer  processes.  The  orga¬ 
nization  and  sequential  implementation  of  these  processes  within 
each  node  is  shown  in  Fig.  3.  These  processes  cumulatively  repre¬ 
sent  the  fundamental  conservations  of  mass,  momentum,  energy, 
and  charge.  The  division  into  the  separate  processes  allows  the 
solution  of  the  ensemble  effect  of  a  number  of  interacting  physics 
via  consideration  of  one  physical  operator  at  a  time,  with  iteration 
allowing  for  coordinated  solution  amongst  the  operators.  This  is 
achieved  by  evaluation  of  the  thermodynamics  process,  which 


Fig.  2.  Sequential  simplification  from  SOFC  stack  to  model  channel-based  domain. 
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Fig.  3.  SOFC  subsystem  process  flow  and  convergence. 


calculates  current  density  as  its  output,  as  the  first  and  last  step  of 
the  operator  sequence.  Convergence  was  then  dictated  by  a  match 
between  the  current  density  in  the  first  and  second  call  to  the 
thermodynamics  process. 

The  thermodynamics  process  utilizes  standard  expressions  of 
voltage  related  to  current  density  and  electrochemical  losses 
modeled  after  the  work  of  Costamagna  et  al.  [34].  The  output  of  this 


process  was  the  current  density  based  on  a  set  of  given  conditions 
within  the  node  and  the  voltage  passed  to  the  process  from  the 
system  controller.  Current  density  was  calculated  implicitly 
through  evaluation  of  the  cell  voltage  via  the  modeled  over¬ 
potential  terms: 


^1  —  P  Activation,  A  V  Activation,  C  P Concentration,  A 

—  P Concentration, C  —  Pohmic  0  ) 

where  E\  is  the  cell  voltage,  Fn  is  the  Nernst  potential,  and  the 
various  7]  represent  the  anode  and  cathode  activation  and 
concentration  overpotentials  and  the  cell  overall  ohmic  over¬ 
potential.  Since  hydrogen  was  assumed  to  be  the  only  electro- 
chemically  active  fuel  in  the  cell,  the  reactions  in  each  half-cell  are 


represented  by 

H2  +02~^H20  +  2e-  (anode  half  -  cell)  (2) 

io2  +  2e~  — >02_  (cathode  half  -  cell)  (3) 

and  the  associated  change  in  Gibbs’  Free  Energy  for  the  reaction 
defined  as 

Agrxn  =  A/lrxn  —  IVriASrxn  (4) 

where 

A/lrxn  =  ^H2o(^Tri)  _  2^°2^Tri^  _  ^H2^Tri)  (5) 


ASrxn  =  ^H2o(^Tri)  ~  ^Oi^Tri)  ~  sH2(JTr[) 

The  Nernst  potential  can  then  be  expressed  as 


En  =  E0 
where 


nF 


hn 


Ph2o 


(6) 

(7) 


E0 


-Agrxn 

nF 


(8) 


where  TTri  is  tri-layer  temperature,  pz  is  the  partial  pressure  of 
component  i,  n  is  the  number  of  electrons  participating  in  the 
reaction  (in  this  case  2),  and  F  is  Faraday’s  constant. 

The  definitions  for  the  loss  terms  themselves  followed  the 
formulation  presented  in  Costamagna  et  al.  with  two  notable 
exceptions.  In  the  current  model,  the  activation  polarization  is 
defined  separately  for  each  electrode  according  to  an  arcsinh 
approximation: 


^Activation ,f  -  anf'  Slnh  '  (2jQ  |  j  (9) 

where  j  is  current  density  and  a  is  the  reaction’s  symmetry  factor, 
assumed  to  be  0.5.  The  calculation  of  the  exchange  current  density, 
jo, i,  was  the  same  as  the  previous  work.  In  addition,  the  definition  of 
the  area-specific  resistance  (ASR)  was  simplified  from  the  previous 
work,  which  included  multiple  factors  for  the  unique  shape  of  the 
particular  SOFC  unit  modeled.  In  the  current  work,  ASR  was  simply 
defined  as 


ASR  —  pAdA  +  Pcdc  +  PElec^Elec  +  PCC^CC,  A  +  PCC^CC,C  (10) 

where  p  denotes  resistivity  of  the  various  materials  and  d  repre¬ 
sents  their  thicknesses.  The  values  of  the  resistivities  of  all 
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Table  1 

Arrhenius  equation  constants  for  modeled  reformation  reactions. 


Reaction 

A 

£a(kJ-mol  ’) 

CH4  +  H20  <->  CO  +  3H2 

6.043el5 

240.1 

CO  +  H20  <-►  C02  +  h2 

1.28e6 

67.13 

CH4  +  2H20  <->  C02  +  4H2 

1.459el5 

243.9 

materials  except  for  the  cathode  were  modeled  as  temperature- 
dependent,  again  following  the  values  formulated  in  [34].  The 
current  work  assumes  an  anode-supported  cell;  thus,  the  thick¬ 
nesses  modeled  in  this  work  were  as  follows: 

dA  =  210  pm,  dc  =  30  pm,  dElec  -  10  pm, 
dec, a  =  125  pm,  dCc,c  =  125  pm 

The  conservation  of  momentum  process  calculated  the  pressure 
drop  in  each  flow  channel  via 

(11) 

where  A P  is  the  pressure  drop,  p  is  the  mixture  density,  V  is  the  flow 
velocity,  DH  is  the  channel  hydraulic  diameter,  and  Ax  is  the  node 
length.  Gas  channel  height  and  width  were  assumed  to  be  1  mm 
and  3  mm,  respectively.  The  flow  was  assumed  to  be  laminar;  thus, 
the  friction  factor  was 


to  the  anode  flow  channel,  thereby  determining  the  change  in 
species  concentrations  and  total  molar  and  mass  flowrate  in  both 
channels.  In  the  current  model,  it  is  assumed  that  hydrogen  is  the 
only  fuel  species  participating  in  electrochemical  oxidation  on  the 
anode  side;  the  effect  of  this  assumption  will  be  discussed  further 
below.  Under  this  assumption,  the  reactions  occurring  in  the  two 
flow  control  volumes  are  as  specified  in  Eqs.  (2)  and  (3).  Extension 
of  Faraday’s  Law  to  a  balance  of  molar  flow  within  each  channel 
then  dictates  that  the  output  mole  fractions  of  all  participating 
species,  after  electrochemistry,  are  described  by 


yH2,  Anode  = 

Nlnlet,AyH2,AJnlet  N[on 

^Inlet,A 

(14) 

yH20, Anode 

_  Nlnlet,AyH20,A,Inlet  —  Njon 

^Inlet,A 

(15) 

yo2, Cathode 

_  ^C,Inletyo2,C, Inlet  —  0-5'Nion 

Nc 

(16) 

yN2, Cathode 

_  Nc,InletyN2,C, Inlet 

Nr 

(17) 

where  N  represents  the  molar  flowrate,  and  Nion  represents 
a  source  on  the  anode  and  a  sink  on  the  cathode  side,  determined 
via: 


where  Re  is  the  flow  Reynolds  number. 

The  reaction  kinetics  process  determined  changes  in  species 
concentrations  due  to  the  action  of  reformation  reactions  and 
water-gas  shift  chemistry.  In  the  current  model,  only  three  reac¬ 
tions  were  assumed  to  take  place.  The  reactions  and  their  rate 
parameters  according  to  the  Arrhenius  equation  are  described  in 
Table  1.  These  model  reactions  are  not  the  elementary  reactions 
expected  to  occur  within  the  cell;  rather,  they  are  summary  reac¬ 
tions  that  represent  the  global  effect  of  a  much  larger  number  of 
elementary  homogeneous  and  heterogeneous  reactions  taking 
place  in  the  channel  and  interacting  with  the  nickel  in  the  anode, 
respectively.  Rate  data  was  obtained  from  a  pilot  investigation 
utilizing  Chemkin,  which  resulted  in  values  slightly  altered  from 
those  reported  in  [6],  after  which  the  reformation  reactions  were 
modeled.  In  addition,  all  reaction  rates  were  inversely  scaled  by  an 
adsorption  rate  constant  described  by  Eq.  (13),  withy*  representing 
the  mole  fraction  of  species  x  and  the  values  listed  in  Table  2. 

DEN  =  ^1  +  /<ads,CO -yco  +  fcads,H2  +  ^<ads,CH4  ^014 

+  kads,H20'^j^  (13) 

The  overall  effect  of  the  DEN  term  is  simply  to  provide  a  mechanism 
of  accounting  for  the  slower  adsorption  kinetics  on  the  Ni  surface  as 
compared  to  the  kinetics  of  the  reformation  reactions. 

The  electrochemistry  and  conservation  of  mass  process  simu¬ 
lates  the  transfer  of  the  oxygen  ions  from  the  cathode  flow  channel 


Table  2 

Adsorption  rate  constants  utilized  in  Eq.  (13). 


Reaction  rate 

A 

Ea  (kj-mol  1) 

^ads,CO 

6.12e— 9 

-82.9 

^ads,H2 

6.65e4 

-38.28 

kads,CH4 

1.77e5 

88.68 

kads,H20 

8.230e-5 

-70.65 

Nin 


2 F 


(18) 


Nc  =  NInlet)C  -  0.5-Nion  (19) 

Finally,  the  heat  transfer  process  determines  the  temperatures  of  all 
four  control  volumes  within  the  SOFC  node.  Since  four  tempera¬ 
tures  must  be  determined,  and  all  are  dependent  on  each  other,  the 
solution  method  for  this  process  was  a  guess-and-check  method, 
similar  to  the  overall  SOFC  subsystem.  Convergence  was  deter¬ 
mined  to  occur  when  the  temperatures  in  all  four  control  volumes 
were  such  that  energy  conservation  could  be  satisfied.  Three 
general  categories  of  thermal  energy  exchanges  were  accounted  for 
in  this  model:  convective  heat  transfers,  conductive  heat  transfers, 
and  thermal  energy  sources  and  sinks. 

Convective  transfers  included  transfers  between  the  anode  and 
cathode  flow  channels  and  both  the  tri-layer  and  current  conductor 
control  volumes.  In  addition,  mass  flows  that  advected  thermal 
energy  through  the  system  were  accounted  for;  these  included 
thermal  energy  entering  and  exiting  the  anode  and  cathode  control 
volumes  through  gas  flow  and  the  thermal  energy  carried  by  the 
oxygen  ion  modeled  as  leaving  the  cathode  gas  flow,  entering  the 
tri-layer,  experiencing  a  near-instantaneous  change  in  temperature 
to  the  tri-layer  temperature,  and  finally  exiting  the  tri-layer  to  enter 
the  anode  gas.  According  to  the  perfectly  stirred  reactor  assump¬ 
tions,  all  thermal  energy  exiting  the  flow  control  volumes  was 
modeled  at  the  converged  temperature  of  the  control  volume 
within  the  node.  For  the  convective  transfers  between  gaseous  and 
solid  control  volumes,  the  convective  heat  transfer  coefficient  was 
modeled  as 


^  = 


Nuj-ki 


D, 


Hj 


(20) 


where  Nuj  is  the  Nusselt  number  correlation  given  the  non-circular 
geometry  of  the  channel  as  well  as  the  aspect  ratio  of  the  channel, 
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Dn,i  is  the  hydraulic  diameter  of  gas  channel  j,  and  /q  is  the  thermal 
conductivity  of  gas  phase  i.  The  correlation  for  Nut  used  in  the 
simulation  was  based  on  empirical  data  provided  in  [35]. 

Conductive  heat  transfers  were  only  modeled  in  the  solid  phases 
between  adjacent  nodes  of  the  same  control  volume  and  also 
between  the  current  collector  and  tri-layer  control  volumes  within 
a  given  node.  In  general,  the  conduction  between  adjacent  nodes  of 
the  same  control  volume  was  modeled  as 


^Condi,i-l  — 


^X-Sec(^-^i-l)  Ai.  r 

dx - dt’£| 


_l<AX-Sec(Ti~TM)  ^ 

Condi, i+l —  ^  ut 

(21) 


where  i  is  the  index  of  the  current  node,  k  is  the  thermal  conduc¬ 
tivity  of  the  given  material,  AX-sec  is  the  cross-sectional  area 
through  which  the  energy  is  transferred,  and  T  is  the  temperature 
at  the  indicated  node.  Conduction  between  the  current  conductor 
and  the  tri-layer  in  a  given  node  was  modeled  as 


EcondiMC  -  WH++n-rcc).dt  (22) 

where  /<Avg,H  is  the  harmonic  average  of  the  two  materials’  thermal 
conductivities,  Ac c-t  is  the  area  of  contact  between  the  current 
collector  and  tri-layer,  T  is  the  temperature,  and  dy  is  the  distance 
from  the  center  of  the  current  collector  to  the  center  of  the  tri-layer, 
in  the  direction  of  travel  of  the  heat  exchange  (e.g.  through  the 
thickness  of  these  two  phases). 

All  thermal  sources  in  the  model  originated  from  the  action  of 
electrochemistry  within  the  cell.  All  of  these  thermal  sources  were 
modeled  to  exist  within  the  tri-layer,  as  they  originate  from  the 
action  of  the  electrolyte  and  electrodes.  Therefore,  the  calculation 
for  all  the  electrochemical  losses’  energy  contents,  £Act,A^Act, 
C,£cnc,A>£cnc,c,  an^  ^ohm  were  derived  from  their  analogs 
presented  earlier.  The  electrical  energy  leaving  the  node  to  the 
external  circuit  was  the  only  sink  and  taken  simply  as 


Eniec  =  IV  dt  (23) 

where  I  is  the  current  and  V  is  the  voltage.  The  final  energy  source 
was  the  entropic  loss  experienced  at  the  standard  potential, 
calculated  as 


^Entropy  —  ^rxn^ion^Tri  ’ 


(24) 


2.3.  Gas  turbine  subsystem  model 

The  gas  turbine  subsystem  utilizes  compressor  and  turbine 
performance  maps  to  determine  the  turbomachinery’s  state  at  each 
timestep.  Convergence  across  this  subsystem  is  determined  by 
satisfaction  of  mass  balance  across  the  gas  turbine,  with  the  map’s 
predicted  mass  flowrate  required  to  match  that  exiting  the  auxil¬ 
iary  combustor.  The  compressor  and  turbine  maps  themselves  were 
modified  versions  of  generalized  performance  maps  readily  avail¬ 
able  from  the  freeware  GasTurb  11.  The  compressor  and  turbine 
maps  were  scaled  to  values  that  represent  the  reported  design 
point  performance  of  the  Solar  Turbines  Saturn  20  gas  turbine 
generator  [36].  This  particular  turbine  was  chosen  based  on 
a  preliminary  study  that  determined  its  power  (1.2  MW)  and  size 
(9.02  m2,  21.74  m3)  were  well-matched  to  the  requirements  of  the 
hybrid  system  needed  in  this  application;  this  development  is 
presented  in  more  detail  in  Part  II  of  this  work. 

Interpolation  was  based  on  the  shaft  speed  at  the  given  timestep 
and  guessed  pressure  ratio.  With  these  two  metrics  known,  it  was 
possible  to  then  use  the  maps  to  determine  the  isentropic  efficiency 
and  mass  flowrate  of  the  compressor  and  turbine.  Map  data  in 


GasTurb  is  provided  along  “(3-lines,”  which  give  sets  of  mass  flow, 
pressure  ratio,  and  efficiency  data  at  various  shaft  speeds  along 
specified  linear  trends  in  the  maps.  Fig.  4  gives  an  example  for  the 
compressor  map;  solid  lines  are  (3-lines  and  dashed  lines  are  lines  of 
constant  shaft  speed.  Similar  relationships  exist  between  the  other 
performance  variables  and  for  the  turbine  performance  as  well. 
Determination  of  state  was  then  a  matter  of  locating  the  p-lines  that 
bound  the  current  state  based  on  pressure  ratio  and  shaft  speed  and 
then  interpolating  between  these  lines  for  mass  flowrate  and 
efficiency. 

The  only  portion  of  the  turbomachinery  performance  that  was 
not  determined  from  the  maps  was  the  shaft  speed.  This  was 
determined  by  a  torque  balance  on  the  shaft,  which  was  itself 
controlled  by  the  load  demanded  of  the  turbomachinery.  Deter¬ 
mination  of  the  shaft  speed  followed 

(25) 

where  Pj,  Pc,  PLoad,  and  Pl0Ss  are  the  power  generated  by  the  turbine 
and  the  power  consumed  by  the  compressor,  load,  and  losses, 
respectively.  Jshaft  is  the  shaft’s  rotational  inertia.  The  loss  term  was 
consistently  set  to  5%  of  the  net  difference  between  turbine  and 
compressor  power.  This  same  relationship  was  used  in  the  SOFC 
average  temperature  control’s  feed-forward  portion  to  predict  the 
load  change  required  for  a  desired  change  in  shaft  speed. 

Finally,  the  auxiliary  combustor  was  also  treated  as  part  of  the 
gas  turbine  subsystem.  The  auxiliary  combustor  accepted  as  input 
the  depleted  anode  fuel  stream,  the  portion  of  the  cathode  stream 
not  recirculated,  and  a  control  signal  for  additional  fuel  flow.  It  was 
assumed  that  complete  oxidation  occurred  within  the  combustor. 
Outlet  species  concentrations  and  temperatures  were  then  calcu¬ 
lated  based  on  these  parameters. 

2.4.  Air  preheater  subsystem  model 

The  structure  of  the  Air  preheater  subsystem  was  very  similar  to 
the  SOFC;  the  heat  exchanger  was  assumed  to  be  a  plate-and-fin 
type  in  counter-flow  and  only  a  single,  centrally-located  channel 
in  the  plate-and-fin  array  was  modeled,  with  spatial  resolution  in 
the  flow  direction.  As  with  the  fuel  cell,  energy  conservation  was 
maintained  in  each  individual  control  volume  (fresh  air,  exhaust 
products,  and  solid)  in  order  to  determine  convergence  of  the 
solution  of  all  three  temperatures  simultaneously.  Determination  of 


l 
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Fig.  4.  Sample  compressor  map;  solid  lines  are  “(3-lines”,  dotted  lines  are  at  constant 
shaft  speed. 
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the  convective  heat  transfer  between  the  solid  (S)  and  a  gas  flow,  i, 
followed  [35],  namely 

^Conv, S->i  =  ^0,i '  (^Plate  +  ^Fin)  ,??Fin '  (^Solid  —  ^Gas)  (26) 

where  ij¥in  is  the  fin  effectiveness,  defined  as 


volume.  These  balances  were  then  solved  iteratively  for  both 
channels,  after  which  the  balance  for  the  solid  phase  was  evaluated 
algebraically. 

2.5.  Control  methodology 


^Fin  -  1  (l  %in)  -A  Fa 

v  '  ''Fin  +  ^Plate 

(27) 

where 

_  tanh  (mrLFin) 

'/Fin  r 

(28) 

/  ho,i-P 

mi  =  \  1 - ~k  - 

V  /<:Sol ' 'VC-Sec, .Fin 

(29) 

^Fin  =  ^Chan  +  W^Fin 

(30) 

where  P  is  the  perimeter  of  the  fin’s  cross  section,  /<soi  is  the 
conductive  heat  transfer  coefficient  for  the  fin’s  material,  Ax_Sec  Fin 
is  the  cross-sectional  area  of  the  fin,  Hchan  is  the  height  of  the  gas 
flow  channel,  and  W¥[n  is  the  thickness  of  the  fin. 

Given  that  the  heat  exchanger  was  modeled  as  counter-flow, 
and  heat  transfer  was  the  only  physics  of  interest  within  this 
subsystem,  a  faster  solution  scheme  was  utilized  than  the  SOFC 
guess-and-check  method.  The  effectiveness-NTU  method  was 
instead  implemented.  An  overall  heat  transfer  coefficient  was 
defined  as 

UA  =  - T - 1 - T -  (31) 

^Con  v,S->C  £conv,S->H 

where  FConv, s^c  and  FConv, s-+h  are  the  convective  heat  transfers  on 
the  cold  and  hot  sides,  respectively,  following  Eq.  (26).  The  net 
transfer  units  (NTU)  was  then  defined  as 

I TA 

NTU  =  — -  (32) 

6HR,Min 

where  CHR  Min  is  the  smaller  heat  rate  from  between  the  hot  and 
cold  sides,  defined  as 

Chr  =  ECp.iCO-JV;  (33) 

1=1 

In  addition,  the  ratio  of  the  heat  rates  was  defined  as 

p  Q-IR,Min  i"3/n 

kHR  =  r - 

<-FIR,Max 

The  effectiveness  of  heat  transfer  within  a  node  could  then  be 
calculated  from  the  above  parameters  according  to 


_  1  -  exp(  -  NTU-(1  -  PHr)) 

1  -  ^HR*exp(  -  NTU(1  -  PHr))  1  } 

where  all  variables  are  as  previously  defined.  Utilizing  this  effec¬ 
tiveness,  the  energy  balance  around  an  individual  node’s  cold-side 
control  volume  would  then  be 

AEGas,C  =  G  •  (Chr, Min'  (^Inlet,H~^Inlet,c))  ~ ^Outlet, C+^Inlet,C  (36) 

A  similar  balance  could  also  be  written  for  the  hot  side,  but  where 
the  convective  heat  transfer  is  out  of,  instead  of  into,  the  control 


The  control  unit  for  the  SOFC-GT  hybrid  system  was  designed  to 
prioritize  safe  operation  of  the  SOFC  before  attempting  to  satisfy 
system  power  requirements.  In  order  to  accomplish  this,  a  cascade 
controller  was  implemented,  with  the  structure  shown  in  Fig.  5.  In 
the  SOFC,  the  major  safety  concern  regards  thermal  stresses, 
induced  by  large  temperature  gradients  and/or  high  operating 
temperature,  which  can  cause  fracture  and  potentially  leaks  of  fuel 
gas.  Thus,  the  two  highest  levels  in  the  cascade  ensure  safe  oper¬ 
ation  is  achieved  first,  with  average  temperature  controlled  by  the 
load  on  the  turbine  shaft  (which  will  affect  the  system  air  flow) 
and  temperature  rise  controlled  by  the  relative  amount  of  depleted 
air  recirculated  on  the  cathode  side.  Efficient  SOFC  operation  was 
then  moderated  by  controlling  anode  utilization  via  manipulation 
of  the  SOFC  voltage.  Finally,  at  the  lowest  priority,  efficient  turbine 
operation  and  satisfaction  of  system  power  demand  were  imple¬ 
mented.  The  former  was  achieved  by  controlling  turbine  inlet 
temperature  through  manipulation  of  the  auxiliary  fuel  flow;  the 
latter  was  achieved  by  manipulating  the  main  system  fuel  flow, 
which  occurs  at  the  SOFC  anode.  Prioritization  was  implemented 
by  setting  bands  in  the  controlled  variable’s  state  of  “safe”  or 
“preferred”  operation,  outside  of  which  control  was  not  allowed  to 
proceed  to  a  lower  level.  When  necessary,  an  additional  target 
setpoint  value  was  specified.  These  specifications  are  provided  in 
Table  3. 

The  majority  of  the  control  was  achieved  through  proportional 
feedback  control.  The  system  net  power  was  the  only  parameter 
observed  to  have  a  steady-state  error  under  P-type  only  control; 
therefore,  integral  control  was  implemented  in  this  instance.  In 
addition,  the  SOFC  average  temperature  control  utilized  a  feed¬ 
forward  signal  to  achieve  greater  stability.  The  physical  parameter 
with  the  largest  direct  effect  on  the  SOFC  temperature  is  the 
flowrate  of  air  through  the  cathode,  which  is  dictated  by  the 
performance  of  the  turbomachinery.  In  a  physical  system  with 
variable  shaft  speed,  air  flowrate  can  be  controlled  by  the  appli¬ 
cation  of  load  on  the  shaft,  as  indicated  in  the  torque  balance  of  Eq. 
(25).  Thus,  in  order  to  implement  the  predictive  capability  required 
by  feed-forward  control,  the  change  in  required  flowrate  in  the 
cathode  would  need  to  be  translated  into  a  change  in  load  on  the 
shaft. 

This  calculation  process  is  embodied  in  Eqs.  (37)— (39),  followed 
by  inversion  of  Eq.  (25)  to  solve  for  load  as  a  function  of  shaft 
speed.  In  order  to  predict  the  necessary  air  flow,  the  heat  load 
generated  by  the  cell  must  first  be  calculated,  per  Eq.  (37)  where 
PMax  is  the  maximum  theoretical  power  (based  on  Nernst  poten¬ 
tial)  at  the  converged  current  density.  This  heat  load  is  assumed  to 
be  the  primary  forcing  of  the  temperature  rise  along  the  cell 
length;  thus,  to  maintain  a  desired  temperature  difference,  an 
equal  amount  of  heat  with  an  equal  temperature  rise  would  need 
to  be  absorbed  by  the  air  flow  in  the  cathode  (it  is  implicitly 
assumed  that  the  flow  in  the  anode  provides  significantly  less 
convective  flux).  Eq.  (38)  then  provides  an  approximation  for  the 
air  flowrate  required.  ATSet  is  the  maximum  allowable  temperature 
difference  along  the  SOFC  channel’s  length  (200  I<),  and  Cp  Ave  is 
the  specific  heat  of  air  at  the  converged  average  SOFC  temperature. 
Finally,  from  the  data  supplied  with  the  compressor  map,  an 
approximate  relationship  between  shaft  speed  and  air  flowrate 
was  developed  as  shown  in  Eq.  (39),  derived  from  a  linear 
regression  of  the  average  relationship  across  all  (3-line  data 
provided. 
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Qcen  =  (^Max  -  ^SOFc)'^Cell 


MAir 


OGen 


Cp.Ave 

Wv^r 


Alset 


JVGT  =  503.3  MAir  + 9880.2 


3.  System  modeling  results  and  discussion 

3.1.  SOFC  subsystem  validation 

f  38) 

^  '  Prior  to  integrating  the  SOFC  Subsystem  with  the  remainder  of 

the  model  code,  this  module  was  validated  against  readily  available 
data  from  a  1996  IEA  project  report  [37].  In  this  report,  several 
modeling  efforts  were  discussed  and  various  SOFC  operating 
(39)  parameters  provided  for  humidified  hydrogen  (90%  FI2;  10%  H20) 
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Table  3 

Controller  targets  and  permissible  operation  bands. 


Controlled  parameter 

Minimum 

value 

Maximum 

value 

Setpoint 

value 

Average  SOFC  temperature 

1067  K 

1179  K 

1123  K 

SOFC  temperature  rise 

-/- 

200  K 

-/- 

Anode  fuel  utilization 

0.7 

0.9 

0.85 

Combustor/turbine  inlet 
temperature 

1223  K 

1523  K 

-/- 

Total  system  net  power 

-/- 

-/- 

Instantaneous  system 
demand 

Table  4 

Comparison  of  SOFC  subsystem  performance  to  available  benchmark  cases. 


Case  A  (H2) 

Case  B  (natural  gas) 

IEA 

Model 

IEA 

Model 

Voltage  (V) 

0.71 

0.71 

0.64 

0.64 

Current  (A-m-2) 

3000  avg 

4203  avg 

3000  avg 

3712  avg 

Power  (W-rrT2) 

2130 

2984 

1923 

2645 

Tanode  (K) 

1332  out 

1293  avg 
1332  out 

1296  out 

1258  avg 
1295  out 

^cathode  (K) 

1332  out 

1280  avg 
1328  out 

1293  out 

1247  avg 
1293  out 

Ttri  (I<) 

1336  max 

1290  avg 
1331  max 
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and  natural  gas  steam  reformate  (26.26%  H2;  17%  CH4;  49.34%  H20; 
2.94%CO;  4.36%  C02)  as  the  feed  fuels,  with  inlet  streams  at  1173  K 
and  1  atm.  The  SOFC  subsystem  was  therefore  run  with  these  same 
conditions,  holding  voltage  at  the  average  values  provided  in  the 
report  and  attempting  to  use  cathode  flowrate  to  match  the  average 
value  of  key  temperatures  reported,  as  well.  Table  4  compares  the 
values  reported  in  the  IEA  document  to  those  obtained  with  the 
SOFC  subsystem.  As  can  be  seen  in  the  data,  fuel  utilization  and 
temperature  features  are  well-matched  between  the  current  and 
previous  works.  The  only  significant  difference  is  in  the  current  and 
power  densities.  However,  this  is  expected  as  the  electrochemical 
model  utilized  in  this  investigation  is  from  a  more  modern  cell’s 
performance  than  those  simulated  in  the  previous  work.  Fig.  6  also 
provides  a  comparison  between  the  IEA  report’s  temperature  and 
current  density  profiles  and  those  obtained  with  the  SOFC 
subsystem,  for  the  H2  case.  Comparing  the  current  model’s  data  to 
those  from  the  IEA  report,  it  can  be  seen  that  key  features,  including 
local  maxima  and  overall  shapes  of  the  trends  are  in  good  agree¬ 
ment.  Note  that  the  improved  performance  from  the  current 
model’s  electrochemistry  is  noticeable  in  that  the  current  density 
does  not  drop  as  sharply  along  the  length  of  the  fuel  cell  channels, 
in  spite  of  similar  maximum  values  at  the  inlet  of  the  fuel  cell. 
Similar  matches  in  trends  were  obtained  in  the  natural  gas  case. 


Fig.  7.  System  efficiency  during  transient  convergence  to  steady  state  for  ail  fuels. 


3.2.  Steady-state  system  design  point  simulations 


With  the  SOFC  subsystem  verified,  the  performance  of  the  entire 
hybrid  system  was  then  investigated,  starting  with  the  expected 
design  point  performance  at  rated  power.  The  system  was  therefore 
initialized  to  a  state  near  the  desired  deign  point  and  then 
controlled  for  an  approach  to  steady-state.  The  target  operating 
state  for  hydrogen  and  natural  gas  reformate  was  to  achieve 
3.5  MW  total  system  power  while  maintaining  an  85%  anode  fuel 
utilization,  an  average  SOFC  temperature  of  1123  K,  a  temperature 
rise  across  the  SOFC  less  than  200  K,  and  a  combustor  temperature 
between  1223  I<  and  1523  K.  For  the  diesel  reformate  case,  the  fuel 
utilization  constraint  was  reduced  to  80%  because  the  higher 
utilization  was  unsustainable  for  this  fuel.  Controller  gains  were 
tuned  such  that  the  performance  could  be  reached  in  the  steady- 
state  simulation  with  reasonable  transient  characteristics; 
however,  they  were  not  optimized  during  this  investigation.  The 
fuel  compositions  for  hydrogen  and  natural  gas  reformate  have 
previously  been  described;  for  diesel  reformate,  the  composition 
was  21.7%  H2;  0.5%  CH4;  17.8%  H20;  10%  CO;  8.8%  C02;  41%  N2.  This 
composition  is  as  observed  from  diesel  autotermal  pre-reformer 
studies  reported  in  [38],  which  utilized  synthetic  diesel  blends  as 
surrogates  for  commercial  diesel. 

The  development  of  the  steady-state  efficiencies  for  all  three 
cases  is  shown  in  Fig.  7.  In  this  work,  the  system  efficiency  pre¬ 
sented  is  a  fuel-to-electric  efficiency  defined  as 
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Fig.  6.  Comparison  of  SOFC  subsystem  spatial  trends  to  available  benchmark  case. 
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where  Psofg  Piurb,  And  Peiow  are  the  SOFC  power,  net  turbine  power 
(referred  to  as  Load  in  the  figures  presented  below),  and  the 
cathode  recycle  blower  power;  LHVS0Fc+Turb  is  the  combined 
heating  value  from  the  flow  of  fuel  into  the  SOFC  anode  channels 
and  the  auxiliary  burner.  All  three  systems  achieve  efficiency  esti¬ 
mates  greater  than  65%,  with  the  natural  gas  reformate  and  diesel 
reformate  cases  achieving  slightly  higher  efficiencies  than  the 
hydrogen  case.  This  trend  in  system  efficiency  was  caused  by 
a  matching  trend  in  the  steady  operating  point  of  the  turbine 
(82.1%,  82.6%,  and  84.1%  isentropic  efficiency  for  hydrogen,  natural 
gas,  and  diesel  cases,  respectively).  The  difference  in  component 
efficiency  was  caused  by  the  combustor  control  method,  which 
allowed  turbine  inlet  temperature  to  acquire  a  range  of  values 
rather  than  holding  all  cases  to  a  common  value.  Simulation  results 
indicated  temperatures  of  1260  K,  1275  K,  and  1330  K  for  hydrogen, 
natural  gas,  and  diesel  respectively.  This  temperature  variation  was 
caused  by  the  lower  SOFC  fuel  utilization  in  the  diesel  case,  which 
allowed  more  fuel  to  be  burned  in  the  combustor,  which  is  a  greater 
thermal  energy  source  than  the  SOFC.  It  should  be  emphasized  that 


the  fuel  reformate  cases  do  not  consider  the  energy  requirements  of 
fuel  conversion;  the  impact  of  the  associated  energy  requirements 
for  this  step  will  be  discussed  in  Part  2  of  this  work. 

For  the  hydrogen  and  natural  gas  cases,  1000  s  was  enough  to 
come  to  steady  state,  but  diesel  required  up  to  10,000  s  (not 
shown).  Additionally,  the  hydrogen  case  was  the  most  stable,  with 
the  natural  gas  reformate  case  exhibiting  some  small  instabilities, 
possibly  due  to  the  competing  action  of  electrochemical  oxidation 
and  reformation  within  the  anode  as  well  as  longer  times  spent  in 
higher  levels  of  the  control  cascade,  leading  to  longer  settling  times. 
Clearly,  the  diesel  reformate  case  exhibited  the  most  dramatic 
transient  behavior,  requiring  long  times  to  settle  at  a  given  level  of 
the  control  cascade  before  proceeding  to  deeper  levels. 

Fig.  8  shows  the  dynamic  response  to  control  decisions  for  the 
diesel  reformate  case.  Panels  a  through  e  display  the  development 
of  the  control  and  manipulated  variables  for  control  levels  1-4, 
respectively  (panels  d  and  e  correspond  to  Level  4  together). 
Panel  f  displays  the  final  species  concentration  development 
through  the  length  of  the  cell  at  10,000  s  of  run  time.  In  panels  a-d, 
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solid  black  lines  represent  the  controlled  variable  on  the  primary 
axis,  dashed  black  lines  represent  the  manipulated  variable  on  the 
secondary  axis,  solid  gray  lines  indicate  the  target  value  for  the 
controlled  variable,  and  dashed  gray  lines  indicate  the  desired 
bounds  of  the  controlled  variable.  There  are  multiple  instances 
through  the  initial  transient  response  when  the  limits  on  average 
SOFC  temperature  were  violated  and  the  control  methodology 
returned  to  this  high  level  before  addressing  the  other  control 
variables.  These  divergences  into  unfavorable  operation  zones 
correspond  to  the  initial  reaction  of  the  system  when  it  enters 
a  deeper  level  of  the  control  cascade.  For  example,  the  first  violation 
in  SOFC  average  temperature  comes  shortly  after  the  first  drop  in 
SOFC  voltage  to  attempt  to  meet  the  utilization  target.  It  can  then  be 
seen  through  comparison  of  panels  a  and  c  that  the  control  of  the 
anode  utilization  is  momentarily  halted  while  the  load  adjusts  to 
meet  the  average  SOFC  temperature  criterion.  This  is  exhibited  in 
panel  c  as  a  short,  horizontal  plateau  in  SOFC  voltage  between  two 
major  drops  in  that  voltage. 


Not  surprisingly,  during  the  times  when  control  returned  to 
Level  1,  the  anode  utilization  itself  exceeded  the  bounds  of  its  desired 
operation;  thus,  once  control  fell  to  this  level  again,  then  another 
change  is  made  in  SOFC  voltage.  As  a  result,  there  is  a  brief  time  of 
approximately  800  s  where  control  seems  to  oscillate  between  Levels 
1  and  3.  Coincidentally,  Level  2  and  the  turbine  inlet  temperature 
restriction  in  Level  4  are  always  within  their  limits  during  this  time, 
therefore  not  requiring  the  attention  of  the  controller.  Finally,  it  can 
be  seen  that  during  these  initial  800  seconds,  the  total  system  power 
is  varying  significantly  around  its  setpoint;  however,  once  Levels  1 
and  3  were  satisfied,  coming  to  compliance  on  the  total  system  power 
was  both  swift  and  apparently  transparent  to  the  higher  control 
levels.  Similar  data  were  collected  for  the  other  fuel  cases;  however, 
the  data  shown  in  Fig.  8  are  representative  and  exhibit  more  of  the 
interplay  between  cascade  levels  in  the  controller  than  the  data  for 
the  other  fuel  cases. 

The  need  for  the  lower  utilization  target  for  the  diesel  case  can 
be  seen  by  close  inspection  of  Fig.  7  and  panels  a-e  of  Fig.  8.  Note 
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Fig.  9.  System  performance  for  diesel  reformate  case  during  development  of  steady  state. 
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that  the  steady-state  solutions  exhibit  a  small  amount  of  oscillation. 
When  the  diesel  reformate  case  was  originally  run  with  a  fuel 
utilization  target  of  85%,  SOFC  voltage  would  be  pushed  too  far 
below  the  limits  of  operation.  As  a  result,  the  control  settled  on 
a  fuel  utilization  of  82%  but  very  low  system  efficiency  and  an 
inability  to  meet  system  net  power  demand.  It  can  therefore  be 
surmised  that  anode  fuel  utilization  in  this  case  near  80%  is  most 
likely  pushing  the  SOFC  Subsystem  to  the  limits  of  its  operational 
range,  imparting  the  observed  instability  on  the  resulting 
performance. 

Panel  f  of  Fig.  8  provides  insight  into  the  fundamental  cause  of 
this  instability.  All  of  the  methane  is  consumed  very  close  to  the 
front  of  the  SOFC  anode  channel,  with  almost  all  of  the  hydrogen 
consumed  by  the  exit.  For  the  case  of  a  general  fuel  reformate 
stream,  methane,  hydrogen,  and  carbon  monoxide  must  all  be 
considered  fuels.  Carbon  monoxide,  much  like  methane,  undergoes 
a  reaction  process  (involving  the  water-gas  shift  reaction)  within 
the  anode  channel  that  gives  rise  to  the  formation  of  fuel  hydrogen. 
Flowever,  as  shown  in  panel  f,  the  water-gas  shift  reaction 
approaches  equilibrium  by  the  time  the  flow  exits  the  cell,  with 
a  high  concentration  of  CO.  This  large  amount  of  remaining  CO 
therefore  limits  the  achievable  utilization. 

It  would  therefore  be  desirable  to  incorporate  a  reaction  model 
that  accounts  for  water-gas  shift,  hydrogen  electrochemical 
oxidation,  and  carbon  monoxide  electrochemical  oxidation  in  order 
to  more  accurately  represent  the  limits  of  operation  on  the  diesel 
reformate.  Flowever,  as  of  the  time  of  this  publication’s  writing,  the 
authors  have  not  identified  a  suitably  detailed  mechanism  that 
accounts  for  the  intricacies  of  the  competition  between  these 
reactions.  Moreover,  it  is  fairly  standard  practice  in  the  literature  to 
assume  that  CO  only  participates  in  the  water-gas  shift  reaction.  It 
should  be  noted  that  the  natural  gas  reformate  case  did  not  exhibit 
this  problem,  most  likely  due  to  the  much  higher  concentration  of 
methane  and  lower  concentration  of  carbon  monoxide. 

Fig.  9  provides  detailed  data  regarding  the  operation  of  the 
individual  power-producing  components  within  the  system  for  the 
diesel  reformate  case.  Panels  a  through  c  present  the  general 
turbine  subsystem  performance,  with  the  operating  line  during  the 
transient  overlaying  the  performance  maps  for  the  compressor  and 
turbine,  respectively  in  panels  a  and  b  and  the  shaft  speed  displayed 
in  panel  c.  Using  the  compressor  map  as  reference,  the  system  was 
first  put  at  a  state  near  the  efficiency  island  at  the  upper-right 
corner  of  the  map.  Progression  from  this  point  to  the  final  steady 
state  then  followed  a  path  of  oscillation,  indicated  in  panel  c. 
Overall,  the  final  steady-state  operation  was  nearer  to  the  efficiency 
island  of  both  components  than  the  initial  state,  with  a  pressure 
ratio  of  6.1.  Within  the  SOFC,  the  current  density  obtained  a  peak  at 
approximately  40%  of  the  total  travel  length  within  the  channel,  at 
645  mA/cm2.  Panels  e  and  f  show  that  Nernst  voltage  dropped 
along  the  length  of  the  fuel  cell,  in  agreement  with  the  rising 
temperature  in  the  same  direction.  Comparison  to  Fig.  8  shows 
current  density  closely  followed  average  SOFC  temperature  and 
Nernst  voltage  closely  followed  the  anode  utilization. 

Table  5  compares  the  final  state  for  all  three  fuel  cases.  There  is 
not  much  variation  between  the  cases  in  the  operation  of  the 
turbomachinery,  with  component  efficiencies  and  shaft  speed  very 


Table  5 

Comparison  of  system  performance  when  operated  on  all  modeled  fuels. 
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J  Ave 

(mA-cm-2) 

V(V) 

ATsofc 

(K) 

Hydrogen 

82.1 

78.2 

14,243 

610.1 

0.743 

151.2 

Natural  gas  reformate 

82.6 

77.1 

14,007 

662.4 

0.677 

187.3 

Diesel  reformate 

84.1 

77.7 

14,149 

523.9 

0.686 

134.1 

close  to  each  other  for  all  cases.  Most  of  the  variation  appears  to 
actually  occur  in  the  fuel  cell  itself,  and  most  apparently  in  the 
voltage  and  the  current  density.  The  current  density  of  the  diesel 
reformate  case  is  markedly  lower  than  the  other  two  cases,  most 
likely  due  to  previous  observations  regarding  the  inlet  composition 
and  the  fact  that  carbon  monoxide  electrochemical  oxidation  is  not 
accounted  for  in  this  model.  In  addition,  the  voltage  of  the 
hydrogen  fuel  case  is  significantly  higher  than  the  other  two  due  to 
its  much  higher  hydrogen  inlet  concentration.  The  natural  gas 
reformate  case  comes  close  to  the  upper  limit  of  temperature  rise, 
possibly  due  to  its  high  current  density,  which  counteracts  a  large 
amount  of  the  endothermic  cooling  provided  by  the  internal 
reformation  at  the  entrance  of  the  cell. 

3.3.  Demonstrative  system  dynamic  response 

In  order  to  investigate  response  to  dynamically  changing  system 
demand,  all  three  cases  were  initialized  to  their  final  states  as 
shown  in  the  previous  steady-state  results.  This  state  was  then  held 
for  an  additional  1000  s  before  a  500  kW  step  in  additional  power 
demand  was  introduced  to  the  controller  logic.  This  power  demand 
was  maintained  for  3000  s,  followed  by  another  step  increase  of 
500  kW,  another  3000  s  holding  period,  a  1000  kW  step  decrease  in 
power  demand,  and  a  final  3000  s  holding  period.  In  a  manner 
similar  to  Fig.  7,  Fig.  10  presents  the  system  efficiency  of  the  three 
fuel  cases.  As  with  the  steady-state  system  results,  the  diesel 
reformate  case  consistently  exhibited  the  highest  system  efficiency, 
with  natural  gas  reformate  very  nearly  matching  the  diesel 
reformate. 

System  efficiency  closely  tracked  the  step  changes  in  demand 
after  the  settling  of  initial  transients,  with  larger  overshoots 
during  load  sheds  than  load  additions.  It  can  also  be  seen  that  the 
effect  is  most  severe  in  the  hydrogen  case.  The  source  of  this  is 
most  likely  thermal  in  nature.  As  the  SOFC  operates,  its  net 
thermal  impact  is  a  production  of  heat.  Increases  in  power  demand 
would  require  higher  temperatures  throughout  the  system, 
including  the  fuel  cell.  Since  the  fuel  cell  prefers  to  generate  heat, 
the  demand  for  extra  heat  to  support  the  additional  demand  is  in 
accordance  with  the  operational  characteristics  of  the  SOFC. 
However,  when  the  system  needs  to  shed  load,  lower  tempera¬ 
tures  are  required  throughout  the  system;  as  the  SOFC  prefers  to 
generate  heat,  rather  than  consume  it,  this  is  in  odds  with  the  fuel 
cell’s  operation.  It  therefore  takes  a  longer  time  for  the  SOFC  to 
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Fig.  10.  Dynamic  response  of  system  efficiency  for  all  modeled  fuels. 
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Fig.  11.  Dynamic  system  performance  during  step  load  changes  for  (a)  hydrogen,  (b)  natural  gas  reformate,  and  (c)  diesel  reformate. 
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Fig.  12.  Electrochemical  dynamic  response  for  (a)  hydrogen,  (b)  natural  gas  reformate,  and  (c)  diesel  reformate. 
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adjust  and  the  impact  on  the  SOFC  performance  is  greater  since  it 
requires  greater  changes  in  the  manipulated  variables  utilized  in 
control. 

In  addition,  the  impact  can  be  expected  to  be  greater  for  the 
hydrogen  case  since  it  does  not  include  any  endothermic  refor¬ 
mation  steps,  which  act  as  a  thermal  sink.  Moreover,  in  the  steady 
state,  the  diesel  reformate  case  requires  ~2.3  times  the  fuel  molar 
flowrate  compared  to  the  natural  gas  reformate  case  because  the 
assumption  of  ATR  in  the  diesel  reformate  case  introduces 
a  substantial  amount  of  nitrogen  into  the  stream,  which  acts  as 
a  significant  diluent  and  heat  sink. 

Fig.  11  presents  the  developed  system  powers  for  all  three  cases 
during  the  transients.  The  hydrogen-fueled  system  responded  the 
fastest,  with  the  smallest  overshoots.  Flowever,  all  three  cases 
provide  acceptable  performance.  The  increased  fuel  feed  require¬ 
ment  of  the  diesel  reformate  case  is  clearly  apparent  in  Fig.  11,  and  it 
can  be  seen  that  for  all  three  cases,  the  system  is  able  to  quickly 
respond  to  the  change  in  demand  by  proceeding  all  the  way  to  Level 
4  of  the  control  cascade.  Upon  the  load  additions,  all  three  fuel  cases 
can  stay  at  this  cascade  level  while  maintaining  compliance  with 
the  other  control  variables.  However,  upon  the  load  shed,  all  three 
cases  leave  this  control  cascade  level,  as  is  exhibited  by  the  brief 
horizontal  plateaus  in  the  Anode  Flowrate  trends.  The  natural  gas 
case  at  first  has  the  most  trouble  of  the  three  fuels.  However, 
investigation  of  both  Figs.  10  and  11  show  that  hydrogen  and 
natural  gas  are  able  to  come  to  steady-state  well  before  the  end  of 
the  settling  time  for  all  three  load  changes.  By  contrast,  the  diesel 
reformate  case  is  still  approaching  steady  state  at  the  end  of 
simulation  after  the  load  shed  (this  is  most  apparent  in  Fig.  10). 

The  longer  settling  time  for  the  diesel  reformate  case  is  an  effect 
of  the  electrochemistry,  and  again  tied  to  the  difficulties  observed 
in  controlling  the  fuel  utilization  when  operating  on  this  fuel.  Fig.  12 
shows  the  electrochemical  performance  of  all  three  cases,  revealing 
that  the  average  current  density  in  the  diesel  reformate  case  is  still 
settling  out  at  the  end  of  simulation,  while  the  other  two  cases 
exhibit  relative  ease  in  coming  to  steady  state.  Fig.  13  shows  the 
voltage  of  the  cell  for  the  diesel  reformate  case  is  nowhere  near 
steady  state  at  the  end  of  simulation,  and  experiences  a  large 
transient  after  the  load  shed  in  an  apparent  attempt  to  bring  the 
fuel  utilization  well  under  control.  In  fact,  for  all  load  changes,  it  can 
be  seen  that  the  fuel  utilization  in  the  diesel  reformate  case  never 
really  reaches  its  target  value,  and  takes  a  long  time  in  its  approach 
to  the  target,  with  significant  voltage  changes  driving  this  process. 
Although  the  settling  time  is  long  for  the  diesel  reformate  case,  at 


Fig.  13.  Dynamic  control  response  for  voltage  manipulation  attempting  to  satisfy 
utilization  constraints  in  the  diesel  reformate  case. 


no  time  does  the  utilization  exceed  the  bounds  of  its  desired 
operation. 

4.  Summary  and  conclusions 

A  model  of  a  simplified  SOFC-GT  hybrid  system  with  cathode- 
side  recirculation  and  the  capability  to  simulate  dynamic  perfor¬ 
mance  as  well  as  provide  spatially  resolved  data  within  the  SOFC 
channels  has  been  successfully  demonstrated.  The  goal  of  the 
model  is  to  analyze  the  potential  for  a  SOFC-GT  system  to  provide 
the  motive  power  in  a  diesel-power  freight  locomotive.  The  design 
of  the  system  model  has  been  built  around  this  capability,  with 
inclusion  of  subroutines  that  can  incorporate  even  more  detailed 
physics  of  diesel  reformation  and  electrochemical  oxidation  of 
reformation  products.  Analysis  of  the  system  began  by  demon¬ 
strating  that  the  performance  of  the  SOFC  subsystem  matches  well 
with  the  published  performance  of  a  benchmark  model  operating 
on  both  humidified  hydrogen  and  natural  gas  reformate.  Simula¬ 
tion  of  the  system  at  its  design  point  has  shown  promising  results 
for  a  system  operating  on  any  one  of  three  candidate  fuels: 
hydrogen,  natural  gas  reformate,  and  diesel  reformate.  Projected 
system  efficiencies  between  65%  and  70%  are  in  agreement  with 
previous  results  in  the  literature  when  taking  into  account  the  fact 
that  the  model  as  presented  does  not  account  for  the  losses  in  the 
reformation  of  the  hydrocarbon  fuels  themselves. 

It  was  found  that  although  all  three  fuels  provide  satisfactory 
performance  and  are  roughly  comparable  at  rated  operation,  the 
use  of  diesel  fuel  reformate  imparts  longer  transients  during 
dynamic  power  changes.  This  was  shown  to  be  related  to  the 
difficulty  in  accurately  and  quickly  controlling  the  fuel  utilization 
for  this  fuel.  The  model  accounts  for  the  large  amount  of  nitrogen  in 
diesel  reformate  (which  acts  as  a  diluent)  as  well  as  the  large 
amounts  of  carbon  monoxide  provided  in  the  inlet  stream.  This 
may  lead  to  the  need  to  critically  evaluate  the  common  assumption 
that  CO  does  not  participate  electrochemically  in  a  typical  SOFC. 
This  work  indicates  there  may  exist  conditions  within  the  cell  at 
which  all  other  reaction  routes  (hydrogen  electrochemistry,  refor¬ 
mation,  and  water-gas  shift)  are  exhausted  and  CO  electro¬ 
chemical  oxidation  may  become  favored. 

Part  II  of  this  work  directly  addresses  the  application  of  this 
system  to  the  locomotive  platform,  concentrating  on  physical 
constraints  as  well  as  performance  along  a  representative  route  in 
the  South  Coast  Air  Basin  in  California.  Utilization  of  modeled  train 
kinematics  as  well  as  physical  characteristics  of  the  representative 
route  will  allow  for  a  more  complete  and  thorough  understanding 
of  the  dynamic  response  capability  of  the  system  than  the  simple 
representative  dynamics  investigated  in  this  work.  Future  work  will 
include  system  design  alternatives,  including  reformer  thermal 
integration.  In  addition,  coke-related  SOFC  degradation  studies  will 
be  carried  out  with  the  goal  to  better  understand  the  system’s 
ability  to  meet  current  service  lifetime  expectations,  typically  on 
the  order  of  decades. 
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